Dynamically correlated regions and configurational entropy in supercooled liquids 
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When a liquid is cooled below its melting temperature, if crystallization is avoided, it forms a 
glass. This phenomenon, called glass transition, is characterized by a marked increase of viscosity, 
about 14 orders of magnitude, in a narrow temperature interval. The microscopic mechanism behind 
the glass transition is still poorly understood. However, recently, great advances have been made in 
the identification of cooperative rearranging regions, or dynamical heterogeneities, i.e. domains of the 
liquid whose relaxation is highly correlated. The growth of the size of these domains is now believed 
to be the driving mechanism for the increase of the viscosity. Recently a tool to quantify the size 
of these domains has been proposed. We apply this tool to a wide class of materials to investigate 
the correlation between the size of the heterogeneities and their configurational entropy, i.e. the 
number of states accessible to a correlated domain. We find that the relaxation time of a given 
system, apart from a material dependent pre-factor, is a universal function of the configurational 
entropy of a correlated domain. As a consequence, we find that at the glass transition temperature, 
the size of the domains and the configurational entropy per unit volume are anti-correlated, as 
originally predicted by the Adam-Gibbs theory. Finally, we use our data to extract some exponents 
defined in the framework of the Random First Order Theory, a recent quantitative theory of the 
glass transition. 



I. INTRODUCTION 

Following the seminal paper of Adam and Gibbs [1], 
the concept of cooperative rearranging regions (CRR) has 
become ubiquitous in the literature on the glass transi- 
tion. In fact, the cooperative relaxation of such regions 
is proposed by many theories to be the elementary rear- 
rangement mechanism taking place in the liquid close to 
the glass transition temperature Tg, and the increase of 
the size of these regions (or dynamical correlation length 
f ) on decreasing the temperature is proposed to be re- 
sponsible for the dramatic slowing down of the dynamics 
around Tg. 

Roughly speaking, the Adam-Gibbs theory proposes 
the existence of CRR of size ^, whose relaxation time 
t(^) is given by: 



to be able to relax, the number of accessible states must 
be larger than a given threshold, let us say rig. Therefore 
^CRR^O > lii(?T.o), which implies by Eq. (2) the exis- 
tence of a lower cutoff ^ on the size of the CRR given by 
KB\'a{no)/pSc{T). Assuming that the dynamical 
decorrelation in the liquid is dominated by the short- 
est relaxation time (smaller CRR), and substituting the 
previous equation in (1), one is left with the celebrated 
Adam-Gibbs relation 
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i.e. it is activated with a barrier proportional to the num- 
ber of units belonging to the CRR (here d is the space 
dimensionality, Kb is the Boltzmann constant, A and 
are system dependent characteristic energy- and length- 
scale respectively). If we indicate by Sc{T) the configu- 
rational entropy per molecule (i.e. the difference between 
the entropy of the liquid and its crystal) , then the loga- 
rithm of the number of accessible states in a CRR, called 
^CHH' ^ central quantity in the AG theory, turns out to 
be 
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where p is the number density of molecules. For a region 



As Sc{T) is experimentally observed to be a decreasing 
function of temperature, and seems to vanish linearly at 
the Kauzmann temperature Tk, both ^ and t{T) are 
predicted to grow and diverge at Tk by the AG theory. 
Note that -by construction-, according to AG theory, the 
" entropy" of the smallest CRR, that dominate the relax- 
ation, is Kb(Tcrr{T) = Sc(T)pi{TY ~ KbHtIo), i.e. 
it is a temperature independent quantity. 

The AG theory had an enormous impact, and the re- 
lation (3) has been shown to be fairly well compatible 
with experimental data. Still, from the theoretical point 
of view, the AG scenario is not firmly established [2], 
mainly because of both the unnatural scaling of r with 
^ in Eq. (1) (one would naturally expect an exponent 
ip < d related to the shape of the CRR-CRR interface), 
and, more important, the predicted values of ^ which 
turn out to be unreasonably small. 

In the last decade, two major advances have been 
made in understanding the relaxation phenomena that 
are behind the Adam-Gibbs picture. Firstly, Kirkpatrick, 
Thirumalai and Wolynes [3] identified a deep analogy be- 
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tween the behavior of supercooled liquids and that of 
a class of mean field spin glass models. These models 
are characterized by the existence of an exponentially 
large (in the system volume V) number M of metastable 
states at low temperature, that give a finite contribution 
Sc{T) = In A/" to the liquid entropy, to be identified 
with the configurational entropy of AG. The mean field 
scenario was then used as a starting point for a nucleation 
theory of supercooled liquids, the Random First Order 
Theory (RFOT) [2, 4, 5]. Secondly, a method to esti- 
mate the size of the CRR in experiments was proposed 
by Berthier at al. [6, 7]. We will now briefly review these 
results. 

As in the AG theory, in the RFOT theory the liquid 

close to Tg is supposed to be a "mosaic state" made of 
CRR of typical radius ^; one assumes that inside a CRR 
the system behaves almost as a mean field system. Then 
one can show [2] that, for a CRR in a state A, the free 
energy cost for nucleation of any possible state B ^ Am 
a droplet of linear dimension r is 

AF^^{r) = -TS,iT)pr'' + rrr (4) 

The configurational entropy turns out to be the driving 
force for nucleation, while T is a surface tension which 
is assumed to be roughly constant around Tg and the 
exponent 6 < d. Note that the free energy barrier for a 
given state is given only by the surfacx; term. The bulk 
term comes from the fact that the number of possible 
different states is exponentially large in the volume of 
the droplet, i.e. there is an cntropic gain in changing 
state inside the droplet. The typical size ^ of the CRR is 
given by the condition AF,g(^) = 0, because for larger 
sizes nucleation inside a CRR is not avoidable and the 
CRR looses its identity [2]. This gives 

The thermodynamic free energy barrier for nucleation 
of a different state inside a CRR is the maximum of 
AF^g(r), which is found in r = (f)^^C and is given 

by 

AF^^ if) oc pi'' T5e(T) = KbT a^^^ (T) , (6) 

i.e. in the RFOT the thermodynamic barrier for nucle- 
ation is given by the total configurational entropy of a 
CRR of typical size Note that in RFOT, using (5) and 

(6), crp^^(T) ~ Sc{T)~'^ and is expected to diverge at 
Tk, while in AG theory it is a constant by definition. 

One of the most interesting open problems in RFOT 
is the relation between this thermodynamic barrier for 
nucleation of a droplet with the relaxation time of the 
system. Usually it is assumed that r ~ e^^-AB^'^y^BT ^ 
e'^cRR^'^') [5]. More generally one can expect that r ~ 
e('^cHH(^))'' , and the exponents -0 and 6 can be adjusted 
to recover the AG relation (3) without imposing Eq. (1). 



While the exponent ip is difficult to compute analytically, 
estimates of the exponent have been obtained by mean 
of instantonic techniques [8, 9]. 

The main problem in RFOT is that the domain size ^ is 
not directly observable; in fact, the quest for a growing 
length-scale in supercooled liquids was for long unsuc- 
cessful, as static correlation functions such as the struc- 
ture factor do not reveal any sign of long range order set- 
ting in on approaching Tg. A major theoretical advance 
in this direction has been achieved in the last decade by 
identifying a family of dynamic [10, 11] and static [2, 12] 
many-points correlation functions, inspired by the mean 
field models, that define correlation lengths which are 
predicted to increase fast on approaching Tg. 

Still these correlations are not directly accessible in 
experiments; nevertheless experimental evidence for a 
growing number of correlated units involved in the re- 
laxation has up to now been obtained by different tech- 
niques, see [13] for a review. 

Very recently, Berthier et al. [6, 7] proposed a very 
general but still simple and direct method to measure 
a "number of correlated units" Ncorr{T) in glass form- 
ing systems. They were able to relate, by a fluctuation- 
dissipation-like theorem, the four-point correlations in- 
troduced in [10] to an easily accessible response function, 
namely the derivative of a dynamic two-point correlation 
(such as, for example, the intermediate scattering func- 
tion) with respect to an external control parameter such 
as temperature or density. Their result can be formulated 
as follows: 

NcorrAiT) = -^^f {maxxT(t)}' , (7) 

where Xrit) = is the temperature derivative of 

a suitable correlation function, and ACp is the config- 
urational heat capacity per molecule at constant pres- 
sure. In fact Eq. (7) is a lower bound for Ncorr,4, but 
one can show [6, 7, 11] that it gives a very good esti- 
mate of this quantity; all the details of the derivation 
can be found in [7] . Moreover one can simplify the anal- 
ysis by assuming that C{t) has a stretched exponential 
form, C(t) = exp(-(VTa(T))'5(^)). Then one has 

^oorrAT) - ^^^^^^ e2 yiWr) ' 

plus two corrections: one involving ^^^ff^, the other com- 
ing from a shift of the maximum of xt that -for large 
stretching- is not found in t = Ta. Both corrections 
can be shown to be irrelevant for the following analy- 
sis, giving an error of the order of 1% on the value of 
Ncorr,4- Note also that the temperature dependence of 
the stretching parameter f3{T) is weak around Tg, espe- 
cially when plotted against In(ra), see e.g. Fig. 2c in [14]. 
Experimental values of Ncorr,A were reported in [6, 7] 
for some prototypical glass-forming systems and it was 
shown that this quantity indeed increases on approach- 
ing Tg. 
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FIG. 1: Original data for NcorrA, computed by i) direct mea- 
surement of the dielectric spectra using Eq. (7), and ii) fitting 
the measured Ta{T) and using Eq. (8), for four different ma- 
terials. In the case of OTP we compare our data with the 
estimate of Ncorr,4 obtained from the VFT fit reported in 
[19]. A discrepancy is observed at small times, while in the 
interesting region of large relaxation times all the results are 
consistent. 

The aim of this paper is to compute Ncorr,i for a large 
number of glass forming materials and compare it to the 
configurational entropy. We hypothesize that Ncorr,A is 
representative of the size of a CRR, i.e. that Ncorr,A ~ 
p^'^; this identification between dynamically correlated 
units and CRR is, in general, non-trivial, as it is possible 
to give examples where it does not hold [7, 11]. Still, 
in the specific case of glass-forming liquids, it is at least 
partially supported by theoretical arguments [2, 7, 11, 
12]. 

Inspired by the AG and RFOT theories, following 
Eq. (2), we define the logarithm of available states in 
a correlation volume as 

As we discussed above, AG theory predicts <Jcrr to 
independent of temperature, while RFOT predicts that it 
is an increasing function of temperature, and predicts re- 
lations between the relaxation time Tq, (T) and (J^rr (T) . 
The aim of this work is to test these predictions. 

From the experimental point of view, there is an impor- 
tant advantage in working with cr^^^ instead of NcorrA- 
In fact in experiments one usually deals with molecular 
liquids or polymers, where the constituents of the liquid 
are complex units and it is not clear a priori what are the 
relevant degrees of freedom which are related to the glass 
transition. Thus, Ncorr,i and Sc have to be normalized 
to a "number of relevant degrees of freedom" (or beads) 
per unit volume. A procedure to define these beads has 
been proposed by Stevenson and Wolynes [15]; still the 
physical meaning of these objects is unclear. On the con- 
trary, cr^^^ represents the number of states accessible to 



a CRR and is independent of the number of relevant de- 
grees of freedom inside this region. In fact, in Eq. (9) the 
dependence on the beads density is cancelled by taking 
the ratio Sc/ACp. 



II. DATA ANALYSIS 

For each analyzed material, we collected original and 
literature data for the relaxation time Tq (T) , the stretch- 
ing exponent P(T), the configurational entropy Sc(T) 
and heat capacity ACp{T) as functions of temperature 
at constant pressure. From these quantities we compute 
Ncorr,A and (JcRR accordiug to (8) and (9), following the 
procedures already detailed in [7]. Some remarks on the 
data analysis are worth at this point: 

1. Sc and ACp are not directly accessible by experi- 
ments. Following an approximation commonly adopted 
in literature, we estimated these quantities from the ex- 
cess entropy of the supercooled liquid with respect to 
the crystal (we used the fits reported in the literature 
or original data obtained for this work). This approxi- 
mation gave rise to several criticisms: for instance, the 
difference of anharmonicity between the crystal and the 
supercooled liquid could reflect in a vibrational contri- 
bution to the excess entropy, see e.g. [16]. One should 
also consider the contribution of the secondary processes, 
when present [17]. All these discrepancies could induce 
an overestimation of Sc up to 50%. Nevertheless, as the 
excess entropy is believed to be often proportional to Sc, 
we used it as an estimation, as previously done in the 
literature. 

2. The determination of the stretching parameter (3 is 
also a possible source of error. Sometimes it has been 
done in the literature by fitting the whole correlation 
function and subtracting the contribution of additional 
fast processes. These procedures are often model depen- 
dent, and can induce an error up to 10% in the value of 
(3. Indeed, in some papers, the /3 of dielectric structural 
relaxation has been underestimated by using a fitting 
procedure on log — log scale, that emphasized the high 
frequency tail of the loss peak. For these cases, /3 has 
been estimated here by using a fitting procedure mainly 
sensitive to the region of the full width at half maximum. 

3. We fitted literature data for Ta{T) by a Vogel-Fulcher- 

B 

Tamman (VFT) law, TaiT) = roe^-^o , in a temperature 
range around Tg where reliable data are available. The 
VFT parameters are reported in table I. For most of the 
systems, we found, as already observed in the literature 
[18, 19], that only one VFT law was not able to fit the 
data in the whole temperature range, and that only in 
a region close to Tg the Adam-Gibbs relation was sat- 
isfied: in these cases only data belonging to this region 
were reported. In some cases we had not access to the 
raw data but only to the VFT fit reported in the litera- 
ture. Note that the characteristic time Tq obtained from 
the VFT fit is not the T ^ cxj limit of Ta{T) as the fit 
is performed in an interval close to Tg where data are 
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FIG. 2: (Top) Ncorr,4 as a function of log(Tc/ro) and (bottom) 
'^CRR ^ ^ function of log(Tc</ro) for the materials listed in 
table I. In the lower panel the spreading of the curves is, 
at fixed log(TQ,/ro), about a factor of 2, while in the upper 
panel is about a factor of 20. The dashed line is log(ra/ro) = 
{o/oo)"*' + z\n{a/ao) + In^ with A = 0.65, ao = 2.86, z = 
1.075, and V = 0.5. 



available. This is consistent since we are interested in 
the "effective" physical attempt frequency close to Tg. 
The derivative d In Tq/c? In T was obtained from the VFT 
fit when it is reliable (see [18] for details). In the cases 
where VFT gave a poor fit, we directly derived the data 
(if dense enough) or we used a polynomial interpolation. 
4. We do not report here a systematic analysis of ex- 
perimental errors; correspondingly error bars are not re- 
ported in the figures. Errors are, generically, of the order 
of 10% on all quantities, mainly coming from the diffi- 
culty in the determination of /3 and Sc as discussed above. 
For some materials different datasets were available: we 
checked that the analysis using different sources gives 
consistent results within these errors. In tables I,II,III 
we report the list of the investigated materials and the 
sources of the data. 

To check the robustness of our results, we measured the 
dielectric permittivity £(w), related to the dipole-dipole 
correlation, of four different substances: OTP (see also 



[7]), DPVC and PPGE, obtained from Sigma- Aldrich, 
and PDF that was kindly provided by Dr. Marian Paluch, 
University of Katowice, Poland. Frequency was spanned 
in the range 0.1 mIIz-3 GHz using a combination of 
Novocontrol Alpha Analyser (up to 10 MHz) and Agi- 
lent Network Analyser (up to 3 GHz). A fine tempera- 
ture scan gives direct access to the temperature derivative 
Xt(w) = ^^^Ff-^, from which we can determine Ncorr,4 by 



using Fq. (7), where time is replaced by frequency. The 
procedure we followed and the results are identical to the 
ones reported in section HI of [7]. Further experimental 
details will be given in [20]. In figure 1 we report the 
values of Ncorr,i for these four substances determined by 
using Fq. (7) and (8): the close agreement between the 
two estimates is a check of the validity of the approxima- 
tions discussed above. Then, for the rest of our analysis, 
we will use Fq. (8) to estimate A"corr,4- 



III. RESULTS 

In figure 2, top panel, we report the parametric plot 
of Ncorr,i{T) versus Ta{T) for the 21 materials listed in 
table I. Our results closely agree with the ones already 
reported in [6, 7]: we observe an increase of Ncorr.i on 
lowering the temperature, consistently with the predic- 
tion of a growing cooperativity in supercooled liquids. 
Note that we added to the plot some materials that were 
not analyzed in [7] ; this is because the liquids studied in 
[7] all have similar values of Sc{Tg) and of VFT parame- 
ters, which then give similar values of Ncorr,i around Tg, 
see table I. We added to the plot materials like selenium, 
TPP and TNB whose configurational entropy at Tg is 
very different from the one of glycerol and other molec- 
ular glass formers. For this reason the spreading of data 
in our figure is more marked than in [7] , about a decade 
at fixed log(Ta/To). Some of these systems are expected 
to severely test any kind of correlation between thermo- 
dynamic and dynamic properties, as they are also known 
as "bad actors" concerning the link between kinetic and 
themodynamic fragility [21]. 

In the bottom panel of figure 2 we show that the data 
collapse much better on a universal curve when (Tcrr is 
plotted instead of Ncorr,i- This result suggests that Tq (T) 
is given by a material-dependent characteristic attempt 
rate tq times a universal function of the thermodynamic 
barrier o'ci?R(^)' i-^- 



log[T„(T)/ro]-.F[a,,,(T)] 



(10) 



with J-[a] an almost universal function. 

To give more robustness to this result one would like 
to add much more materials to the plot. Unfortunately, 
for most materials calorimetric data and/or the stretch- 
ing exponent are reported only at Tg, and the previous 
analysis is not possible. We thus resort to a weaker 
test: for most materials tq ^ 10""'^^ s, therefore the 
glass transition temperature Tg is roughly defined by 



5 



10 



10 



\9 

o. , 



• materials of table 1 
(D moleculars (table 2) 
o polymers (table 3) 
--FIT: y=310/x 



X 



10 s (T ) 



10 




FIG. 3: Plot of NcorrA(Tg) as a function of Sc{Tg) (in units 
of R) for the 45 materials listed in tables I,II,III. The dashed 
line is a fit to the inverse correlation predicted by Eq. (11). 
In the inset the same plot is reported in linear scale. 



\og[Ta{Tg)/To] = 17. If Eq. (10) holds, we expect 



^CRR C^g 



Kb 



NcorrA{Tg) ^ , (11) 



which implies an inverse correlation between Sc and 



Nr.. 



at Tg. These quantities can be easily computed 



from the values of Sc{Tg), ACp{Tg), (3{Tg) and fragihty 
m = d^og^ \T=Tg ; that are reported in the literature 
for many more materials, see tables II and III for their 
list and the references. We can then test the prediction 
(11) on a larger set of 45 materials. In figure 3 we plot 
Ncorr,4{Tg) as a function of Sc(Tg); the plot is compatible 
with an inverse correlation of these two quantities. 

Note that from Eqs. (11), (8) it follows easily that 
Sc{Tg)f3'^m'^ / ACp{Tg) = const. Using the well-known re- 
lation m/17 ~ ACp{Tg)/Sc{Tg), that follows from the 
Adam-Gibbs relation (3), we obtain p'^m = const. The 
latter relation has already been derived in [5] in a differ- 
ent way and seems to be well verified in real materials, 
which is a nice consistency check of our results. 



IV. EXPONENTS 

Having stated our main experimental results, Eqs. (10) 
and (11), we can turn to a more detailed comparison with 
RFOT. First we note that the experimentally observed 
increase of o'crrC^) lowering the temperature is not 
compatible with AG theory, while it fits well into the 
RFOT scenario. Then we would like to extract the val- 
ues of the RFOT exponents from the experimental data: 
we anticipate that this is a very difficult task given the 
large experimental errors already discussed and the prob- 
able importance of preasymptotic effects. Still we observe 
that, identifying Ncorr,4 oc p£,'^ and using Eqs. (5), (6), 



FIG. 4: Parametric plot of Ncorr,4,iT) as a function of 
'^crrX^) materials in table I. The dashed line has 

slope 4/3 and the dot-dashed line has slope 1. 



n: 



e/d 



corr,4 



(12) 



RFOT predicts the relation 

In figure 4 we plot Ncorr,i as a function of cTc-hh ^'^^ 
materials of table I. The power-law relation predicted by 
RFOT is very well verified, and the resulting values for 
the exponent are in the range 2.2-^2.5, see table I. Note 
that RFOT usually assumes 9 = 1.5 [5], while istantonic 
calculations give 6* = 2 [8, 9]. This result is then quite 
puzzling; note however that the relation Ncorr,A ^ is 
not well established [7] and this could affect the result 
for 6. 

We can also try to estimate the exponent ip that re- 
lates the relaxation time to ct^rr: ~ tob'^crr. To 
account for preasymptotic effects, we choose a specific 
form J^[a] = {ajooY' + zln(cr/cro) + InA [7] in Eq. (10) 
to fit the curves in the bottom panel of figure 2. We 
obtain as a best fit a value -0 ~ 0.4 ^ 0.5; note however 
that, by changing the other parameters, one can obtain 
quite good fits for a range of values of -0 = 0.3 1.5. 



From log(Ta/To) ~ ^c„k ^'^d Eq. (12) one obtains 



log- 



To 



5e(T)- 



(13) 



which is consistent with the Adam-Gibbs relation (3) 
only if = 1. Assuming this and "0 ^ 0.4 0.5, we get 
6* ~ 2 2.15 which is consistent with the values reported 
in table I. The coincidence of these two estimates of B 
might support the robustness of the indication that its 
value should be close to or slightly larger than 2. 

The result for 0; ~ 0.4-^0.5 is less reliable. Still, we ob- 
tain a strong indication that < 1, that implies that the 
relaxation time diverges slower than the exponential of 
the thermodynamic barrier. This suggests the existence 
of relaxation paths which are more efficient than simple 
nucleation of a random state inside a CRR. A theoretical 
description of such processes is still lacking. 
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V. CONCLUSIONS 

We collected original and literature data on a set of 
45 glass-forming materials and, exploiting the methods 
of [6, 7], we show that i) the size of the CRR is inversely 
correlated with the configurational entropy per unit vol- 
ume at T = Tg, Eq. (11) and ii) more generally that 
the relaxation time seems to be a universal function of 
the configurational entropy of a CRR, Eq. (10), in a wide 
range of temperatures around Tg. We compared these re- 
sults with theoretical predictions and found good agree- 
ment with the RFOT scenario. Wo also gave an estimate 
of the exponents 6 and -0 of RFOT and found unexpected 
but still reasonable values of these exponents. 

Wc wish to stress again that in our analysis we made 
the strong hypothesis that A^corr,4) as defined in [7], is 



representative of the size of a CRR, i.e. that Ncorr,^ oc 
p^'^. A different scaling would not only change our esti- 
mates for the RFOT exponents, but change the definition 
of (TcRR ^^'^ ultimately destroy the scaling of the curves 
in the bottom panel of figure 2. The consistency of our 
results seem to support the validity of this assumption, 
see however [7, 11] for a more detailed discussion. 
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TABLE I: Summary of data at Tg and references. When more databases were present, we compared the data and verified the con- 
sistency. [•] indicates original data from this work. Most of the dynamic data come from dielectric relaxation experiments, except 
for Ge02 and Si02 (viscosity was used instead of r), for B2O3 and ZnCb (photon correlation spectroscopy), selenium (mechan- 
ical relaxation). Temperatures are in K, times in s, entropy and specific heat in units of R. Note: PC=Propylene Carbonate, 
TNB=TriNaphthyl-Benzene, OTP=0-Terphenyl, PDE=Phenolphtaleindimethyleteher, DPVC=Diphenyl-vinylidene Carbon- 
ate, PPGE=Polyphenylglycidylether, PG=Propylene Glycol, MTHF=2-Methyltetrahydrofuran, TPP=Triphenylphosphite 



Name 


Ta 


log TO 


B 


To 


Sc 


ACp 


13 


rn 


^corrA 


e 


Ref dyn. 


Ref. calor. 


PC 


155 


-14.8 


467 


128 


1.86 


9.10 


0.7 


94 


341 


2.2 


[22, 23] 


[24, 25] 


TNB 


344 


-18.0 


1620 


264 


6.22 


18.1 


0.56 


86 


92 


2.4 


[26, 27] 


[5, 24, 28-31] 


OTP 


244 


-14.5 


684 


202 


2.67 


13.8 


0.55 


97 


148 


2.3 


[18, 19, 32, 33],[.] 


[19, 25, 34] 


PDE 


294 


-20.7 


1793 


215 


4.22 


13.2 


0.73 


84 


207 


2.4 


[33, 35, 36],[»] 


[•] 


DPVC 


251 


-19.0 


1243 


192 


3.24 


12.3 


0.67 


89 


208 


2.2 


[37],[.] 


[•] 


PPGE 


258 


-15.1 


708 


217 


4.97 


23.0 


0.49 


107 


85 


2.3 


[38],[.] 


[39] 


m-toluidine 


185 


-14.9 


519 


154 


1.77 


11.2 


0.57 


102 


217 


2.2 


[40, 41] 


[42] 


PG 


168 


-12.9 


708 


120 


2.28 


8.06 


0.72 


53 


129 


2.1 


[43] 


[44-47] 


GeOa 


816 


-13.8 


9732 


199 


1.08 


0.75 


1 


21 


413 


2.4 


[48-50] 


[21, 50] 


Si02 


1452 


-13.8 


14562 


530 


0.65 


0.37 


0.7 


25 


579 


2.5 


[48-50] 


[48] 


ZnCl2 


385 


-12.8 


1647 


274 


0.63 


1.89 


0.71 


51 


499 


2.2 


[51] 


[45, 52] 


SBromo-pentane 


108 


-12.9 


374 


83 


2.81 


9.09 


0.62 


64 


124 


2.3 


[19, 53] 


[19] 


MTHF 


90 


-17.3 


406 


69 


2.81 


8.91 


0.62 


83 


217 


2.3 


[19, 53] 


[19] 


n-propanol 


100 


-10.5 


386 


70 


2.43 


6.09 


0.62 


41 


78 


2.3 


[19, 53] 


[19] 


Salol 


221 


-15.9 


823 


175 


3.46 


13.2 


0.58 


86 


135 


2.3 


[19, 53] 


[19] 


Butyronitrile 


95 


-11.9 


326 


72 


0.82 


4.84 


0.77 


56 


285 


1.8 


[54] 


[55] 


TPP 


204 


-17.5 


704 


169 


4.18 


18.2 


0.51 


111 


125 


2.3 


[25, 56] 


[25] 


Selenium 


309 


-15.5 


1077 


248 


0.83 


1.83 


0.42 


88 


535 


2.6 


[57, 58] 


[45, 59] 


Glycerol 


188 


-14.7 


973 


130 


3.18 


10.0 


0.68 


54 


97 


2.1 


[53, 60, 61] 


[24, 25] 


Toluene 


116 


-15.1 


328 


97 


1.49 


8.71 


0.55 


103 


267 


2.2 


[62, 63] 


[25] 


B2O3 


553 


-10.7 


2540 


353 


3.18 


3.25 


0.62 


35 


107 


2.5 


[64, 65] 


[16, 21, 31, 45, 46] 



TABLE II: Data at Tg and references for the molecular liquids used in fij 


jure 3. Same units 


as in table I. 


Name 


T, 


5c 


ACp 


P 


m 


N corr,4 


Ref. dyn. 


Ref. calor. 


Sorbitol 


268 


4.41 


28.8 


0.48 


128 


94 


[53, 66] 


[21] 


Ethylbenzene 


115 


2.82 


9.67 


0.68 


58 


115 


[21] 


[67] 


Isopropyl Benzene 


126 


2.52 


10.2 


0.56 


90 


179 


[41, 53] 


[21, 25] 


Triphenylethene 


248 


5.46 


14.5 


0.5 


91 


102 


[21, 68] 


[16, 21] 


Ethylene Glycol 


153 


2.04 


7.58 


0.78 


52 


156 


[21, 69] 


[16, 21] 


Ethanol 


97 


1.63 


5.46 


0.7 


55 


195 


[21] 


[16, 21] 


3-methylpentane 


77 


3.07 


8.66 


0.62 


56 


100 


[21, 53] 


[16, 21] 


Diethyl phthalate 


180 


3.03 


15.3 


0.64 


78 


117 


[70] 


[16] 


a-phenyl-cresol 


220 


3.34 


15.4 


0.53 


83 


90 


[53, 71] 


[71] 


Glucose 


309 


2.97 


17.3 


0.37 


115 


74 


[21, 72] 


[21, 45, 46, 73] 


Indometacin 


318 


4.55 


19.8 


0.59 


75 


71 


[21, 74] 


[21] 
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TABLE III; Data at Tg and references for the polymeric liquids used in figure 3. Same units as in table I. For some systems 
reported by [17], a particular procedure, introduced there, was used to subtract from the excess entropy the contribution of 
the secondary relaxation. 
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Sc 




13 


m 




Ref. dyn. 


Ref. calor. 


PVC 


354 


0.28 


2.33 


0.27 


191 


818 


[53, 75] 


fl7, 75] 


PET 


342 


1.08 


9.36 


0.48 


156 


430 


[53, 75] 


fl7 751 


a-PMMA 


378 


0.58 


3.61 


0.37 


145 


572 


[53, 75] 


[17 751 


PS 
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0.60 
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0.35 


143 


528 


[53, 75] 


fl7, 75] 

[-1- 1 , \ 


PP 
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0.37 


137 
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[53, 75] 


fl7, 75] 
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[17, 75] 


PEN 


390 


0.85 


10.1 


0.48 


140 


320 


[17, 53] 


[21] 


PEEK 


419 


0.54 


10.4 


0.32 


280 


554 


[17, 53] 


[17] 


Pcarb 


420 


1.22 


7.21 


0.35 


132 


212 


[17, 53] 


[17] 



[56] B. Schiener, A. Loidl, R. V. Chamberlin, and R. Bohmer, 

J. Mol. Liq. 69, 243 (1996) 
[57] R. Bohmer and C.A. Angell, Phys.Rov.B, 48, 5857 (1993) 
[58] CM. Roland, P.G. Santangelo, D. J. Plazek and K. M. 

Bernatz, J. Chem. Phys., Ill, 9337 (1999) 
[59] S. S. Chang and A. B. Bestul, J. Chem. Thermodyn. 6, 

325 (1974) 

[60] N. Menon, K. P. O'Brien, P.K. Dixon, L.Wu, S.R. Nagel, 
B.D. Williams and J.P. Carini, J. Non-Cryst. Sohds, 141, 

61 (1992) 

[61] A. Schonhals, F. Krcmcr, A. Hofmann, E.W. Fischer, E. 

Schlosser, Phys. Rev. Lett., 70, 3459 (1993) 
[62] A. Doss, G. Hinze, B. Schiener, J. Hemberger, R. 

Bohmer, J. Chem. Phys. 107 (1997) 1740. 
[63] K.L. Ngai, S. Capaccioli, Phys. Rev. E 69, 031501 (2004) 
[64] D. L. Sidcbottom, R. Bergman, L. Borjcsson, and L. M. 

Torell, Phys. Rev. Lett. 71, 2260 (1993) 
[65] D. L. Sidebottom, B. V. Rodenburg, and J. R. 

Changstrom, Phys. Rev. B, 75, 132201 (2007) 
[66] S. Hensel-Bielowka, M. Paluch, J. Ziolo, and C. M. 

Roland, J. Phys. Chem. B, 106, 12459 (2002) 



[67] O. Yamamuro, L Tsukushi, A. Lindqvist, S. Takahara, 
M. Ishikawa, and T. Matsuo, J. Phys. Chem. B 102, 1605 

(1998) 

[68] B Jakobsen, K. Niss and N.B. Olsen, J. Chem. Phys., 

123, 234511 (2005) 
[69] S.S.N. Murthy, J. Phys. Chem. B, 101, 60, 43 (1997) 
[70] S. Pawlus, M. Paluch, M. Sekula, K. L. Ngai, S. J. 

Rzoska,l and J. Ziolo, Phys. Rev. E, 68, 021503 (2003) 
[71] S. S. N. Murthy, A. Paikaray, and N. Arya, 102, 8213 

(1995) 

[72] Gangasharan, S.S.N. Murthy, J. Phys. Chem., 99, 12349 
(1995); Gangasharan and S.S.N. Murthy, J. Chem. Phys. 
99, 9865 (1993) 

[73] J.D. Stevenson, P.G. Wolynes, J. Phys. Chem. B, 109, 

15093 (2005) 

[74] L. Carpentier, R. Decressain, S. Desprez, and M. 

Descamps, J. Phys. Chem. B 110, 457 (2006) 
[75] C. M. Roland, P. G. Santangelo, and K. L. Ngai, J. Chem. 

Phys., Ill, 5593 (1999) 



